VBS State Induced by Impurity Frustration in CrSNi 
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We provide a physically meaningful picture of the nature of the ground state of the CrSNi com- 
pound in the regime where it is a spin singlet. According to this picture, the anisotropy of the Ni 
atom in the Cr ring induces a dimerization in the molecule that makes the ground state to stabilize 
in a Valence Bond Solid phase of virtual spins. We characterize rigorously this phase by means of 
a particular non-local order parameter denoted the generalized string order parameter. In the com- 
pletely antiferromagnetic regime, the system becomes frustrated. We have performed a numerical 
real-time evolution study of the correlations between the spin of the Ni impurity and the rest of 
the spins in order to show the reaction of the system under this frustration. We provide results that 
manifestly show the agreement with the Valence Bond Solid picture introduced. 

PACS numbers: 75.10.Jm, 75.50.-y, 75.10.-b, 74.20.Mn 



I. INTRODUCTION 

Recently, the physics of ring-shaped molecular magnets 
with antiferromagnetic interactions and an odd num- 
ber of interacting spin centers (e.g. paramagnetic ions) 
has attracted a great deal of interest since they provide 
emblematic examples of systems where spin frustration 
effects due to quantum magnetism play a major role. 
Moreover, they have been sinthcsizcd and studied exper- 
imentally [U, H, i, i, @. Specifically, we shall concen- 
trate on the anomalous magnetic properties in some het- 
erometallic odd spin rings p, 0, S namely, chromium 
rings. The system comprises eight chromium(III) ions 
with spin S" = | each, and one nickel(II) ion with spin 
S ^ \. The magnetic properties of this first odd-member 
antiferromagnetic ring has been investigated with elec- 
tron paramagnetic resonance (EPR) and its spin frus- 
trated properties has been visualized by means of a 
Mobius strip. In this paper we propose an alternative 
and complementary picture of the ground state of this 
CrSNi ring molecule using valence bond states (VBS) 
[13, [m, [ia [il Q Hi, [11, M [13 of virtual spins which 
are used to represent the spins S = | , 1 of the real con- 
stituent ions. We will show that the particular bond pat- 
tern acquired by these VBS states is a consequence and 
manifestation of the spin frustration in the odd CrSNi 
ring molecule. 

Molecular nanornagnets are fascinating new magnetic 
materials [l^ [H, [2l[. They appear in a large vari- 
ety of compounds with many different properties. We 
shall focus on antiferromagnetic compounds of bimetal- 
lic rings. These molecules are ideal candidates to study 
the physics of simple but nontrivial spin models like the 
AF Heisenberg interaction. The key point here is that, 
these molecules show very interesting finite-size quantum 
many-body effects which are typically overlooked in other 
studies of the Heisenberg model where the main focus is 
to achieve the thermodynamic limit (number of spins go- 
ing to infinity). In those studies, the small finite-size ef- 
fects are considered spurious effects that vanish for larger 
and larger systems, which eventually may show some sort 



of universality, if that is the case. Quite on the contrary, 
the nice thing of these small molecules is that we can vary 
its size and coupling constant strenghts such that the 
finite-size effects become some real property that can be 
addressed experimentally, theoretically and numerically. 
Some intersting examples of these small quantum effects 
that we study in this paper are level crossing, change in 
the nature of their the ground state (e.g., from spin sin- 
glet to spin triplet or higher), existence of excited states 
very close to the ground state, etc. 

In Sect In] we introduce a Heisenberg Hamiltonian to 
describe the interactions between the two types of ions 
in the bimetallic compound CrSNi, where the Ni ion 
plays the role of an impurity within an homogeneous 
chain of Cr ions with the shape of a ring molecule. In 
fig. [T] we present the energy spectrum of this Hamilto- 
nian obtained numerically with an appropriate Lanczos 
technique. In Scct lllH we first provide the VBS-picture 
for the CrSNi ring molecule based on a strong coupling 
limit in the Ni impurity coupling. This is the origin of 
the spin frustration in the system. In order to support 
this VBS-picture, we provide numerical results for a gen- 
eralized string order parameter that is able to detect the 
type of virtual bond structure. In Sect llVI we study the 
frustration effects in this ring-molecule and test the asso- 
ciated VBS state picture by means of a numerical study 
of time-evolved correlation functions of spin-spin oper- 
ators for different ions in the nano- molecule. Sect|V] is 
devoted to conclusions. 



II. HAMILTONIAN DESCRIPTION OF THE 
CRSNI MOLECULE RING 

Magnetic molecules are emblematic instances of an en- 
semble of non-interacting quantum systems embedded 
in a condensed matter environment. The synthesis of 
molecular magnets has undergone rapid progress in re- 
cent years. Each of those quantum systems are identi- 
cal molecular units that can contain as few as two and 
up to several dozens of paramagnetic ions (spins). In 
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our case, they correspond to CrSNi ring molecules. This 
molecule is one of the many relevant molecules contain- 
ing transition-metal ions whose spins are so strongly- 
exchange coupled that when the temperature is low 
enough, their behaviour is like single-domain particles 
with a certain total spin [l^, [IS, H!! • 

Macroscopically, these materials appear as crystals or 
powders. Nonetheless, their intermolecular magnetic in- 
teractions are utterly negligible when compared to their 
intramolecular interactions. Thus, measurements of their 
magnetic properties reflect mainly ensemble properties of 
single molecules. There are two major advantages in the 
research on these molecule aggregates. Firstly, the out- 
standing degree of accuracy by which their magnetic dy- 
namics can usually be modeled. Secondly, the opportu- 
nity to chemically engineer molecules possessing desired 
physical properties [li, [l^, HH . In our case, the interest 
relies on the study of the Heisenberg model in situations 
that are usually discarded when studying that model in 
infinite one-, two-, and three-dimensional systems. Our 
studies will reinforce the idea that such spin arrays yield 
qualitatively new physics caused by the finite size of the 
system. 

The CrSNi compound belongs to a wider family of 
molecular rings. In molecules with a small number of 
ions, there exist big differences in their physics depend- 
ing on each particular compound. A first major difference 
is related to the number of its ions: odd or even. In this 
regard, the implications are mainly twofold: first, having 
a difference of one ion in the same family of molecules, 
can cause the molecule to have a neat magnetic mo- 
ment or not, therefore changing its magnetic properties 
drastically. Second, and interestingly enough, a molecu- 
lar ring with completeley antiferromagnetic interactions 
between nearest neighbors ions can be a candidate to 
present quantum spin frustration properties if the num- 
ber of atoms is odd, while this effect will not generally 
be present for even member rings in a given family com- 
pound. 

As it happens, in the majority of these molecules the 
localized single-particle magnetic moments of the ions 
couple antiferromagnetically. Then, their spectrum is de- 
scribed rather well by the Heisenberg model with very few 
parameters because of the high symmetry of the molec- 
ular configurations. These coupling parameters corre- 
spond to isotropic nearest neighbor interaction some- 
times augmented by anisotropy terms. 

The CrSNi compound is one of the first antiferromag- 
netic odd member rings which has been artificially sinthe- 
sized. The results of its magnetic properties are inter- 
preted within the framework of a spin Hamiltonian ap- 
proach and they nicely fit the pattern of the energy levels 
obtained by inelastic neutron spectroscopy. There exists 
also reports on its magnetic and spin frustration effects 
7|. In view of these properties, it has been proposed 
l| that the behavior of this molecule can be properly ex- 
plained with a nearest neighbors Heisenberg model where 
only two different microscopic couplings play a role: one 
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FIG. 1: Energy spectrum of some of the lowest lying states 
in each sector with well defined total spin. In the coupling 
region shown, the spectrum is highly dependent on the values 
of J and J'. In particular, in the antiferromagnetic region 
J' > K the ground state is a triplet for J' < 1.5 K and a 
singlet elsewhere. Inset: for higher values of J' the ground 
state is a singlet and the first excited state is a triplet very 
close in energy. 



is the coupling that parametrizes the strength of the in- 
teraction existing between the Ni and the neighbor pair 
of Cr. The other one is the coupling that takes into 
account the interaction between the Cr-Cr pairs, which 
can be considered the same for each pair. The easy-axis 
anisotropy term is reported to be very weak how as to 
play any role. Therefore, the Hamiltonian that we shall 
study has the following form: 



i=l 

+ J' [SNi-Scr(f) + SNi-Scr(8)], 



(1) 



where, for convenience, the Cr atoms have been num- 
bered from 1 to 8, being these latter the two neighbors 
of the Ni atom. Notice that, since the spin of the niquel 
is equal to Sni = f and the spin of the Cr atoms is 
Scr = 3/2, the total spin of the molecule must be inte- 
ger. 

In fig. [1] we have plotted some lowest lying energy lev- 
els of this Hamiltonian. In Sect. llVl we shall explain how 
these numerical results have been obtained with a multi- 
target Lanczos method. It can be observed that, in some 
regime of the coupling constants J' and J, the energy lev- 
els are highly braided and, as a consecuence, the ground 
state has different value of the total spin depending on 
the exact value of these couplings. However, as can be 
seen from the inset of this figure, for large values of J' 
the ground state is always a spin singlet with a triplet 
state very close in energy above it. In this work, we will 
restrict our study to the antiferromagnetic (J' > 0) re- 
gion and in particular to the region where the ground 
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state is a singlet. From fig. [T] we see that this area cor- 
responds to J' > 1.5 K, while the region < J' < 1.5K 
is characterized by a ground state with total spin equal 
to one (for convenience the computations in fig. [T] have 
been done with a fixed value of J = 16 K). 

The interest in the domain where the ground state is 
a singlet comes not only from the fact that it spans the 
most extension in the antiferromagnetic area, but also be- 
cause the physics of the real CrSNi seems to be in agree- 
ment with a regime close to J = 16 K and J' 70 K, 
with a non-magnetic ground state. Therefore, the inter- 
est of our study relies on the fact that it can provide new 
insights into the physics of a not so well-known state of 
matter but with a well defined connection with experi- 
ments in real compounds. 



III. IMPURITY INDUCED VBS PICTURE 

A valence bond solid is a particular quantum many- 
body state that can be understood as follows: given a 
system of real particles with total spin 5, we can split 
each one of them into 25 virtual particles of spin S = 
1/2. In order to recover the original spins, we enforce 
these virtual particles to couple (i.e., symmetrize) among 
themselves in order to give the original spin S particles. 
To create now a wavefunction with total spin equal to 
zero, we make singlets (i.e., antisymmetrize) out of every 
pair of virtual particles. 

We will denote each of these singlet pairs between vir- 
tual particles as a bond. There are a lot of different 
possible ways to fix the bonds between all the virtual par- 
ticles and. in general, the total wave-function may have 
contributions from all these configurations. There ex- 
ist however some physical situations in which only some 
particular bond configurations, out of the whole possible 
set, take part in the wave-function: some systems have a 
major contribution coming only from one particular bond 
arrangement. These systems are commonly dubbed bond 
crystals. It may also happen that there exist not one but 
some few bonds configurations whose weights are domi- 
nant in the total wave-function. In this case the system 
is called a resonating valence bond sohd (RVBS). There 
exist yet another kind of more disordered states, denoted 
(m,n)-VBS, which we shall see that describe properly 
the ground state of the CrSNi in the singlet region. A 
general (m, n)-VBS state is built by forming bonds only 
between virtual spins belonging to neighbor real parti- 
cles, with the numbers m, n satisfying m + n ~ 2S and 
S being the spin of the real particle. 

These states are usually translationally invariant (with 
the (to, rt)-VBS notation this means that to = ti) in those 
systems where the Hamiltonian possess this symmetry. 
There also exist dimerized (m, n)-VBS states that have 
been shown to appear in systems where the full transla- 
tional symmetry of the Hamiltonian has been partially 
broken such that still exist an enlarged unit cell. Typ- 
ically, this effect can be obtained by introducing an ex- 




FIG. 2: Valence Bond Solid picture of the CrSNi ring 
molecule. The bonds established between the two virtual 
spins in the Ni and the two virtual spins in the neighbor Cr 
atoms forces a dimerized pattern in the rest of the chain. This 
VBS configuration becomes the dominant one in the strongly 
coupled impurity limit J' ^ J. 



ternal dimerization coupling constant in the Hamiltonian 
that still preserves some periodicity. 

Our main result in this section is that there is another 
mechanism to provide such dimerized (m, n)-VBS states 
in the CrSNi ring molecule and whose success precisely 
resides in the existence of an impurity within an homo- 
geneous system. To understand this mechanism in the 
particular case of the CrSNi, we resort to the strong cou- 
pling limit where the antiferromagnetic interaction be- 
tween the Ni and its neighbors is much larger that the 
interaction among Cr pairs. As shown in fig. [2l the two 
virtual spins comprising the Ni will be likely to form 
bonds with the virtual particles in the neighbor chromi- 
ums to satisfy their antiferromagnetic constraints. The 
rest of the virtual spins left will then tend to form similar 
bonds with their neighbor partners, giving as a result a 
dimerized non translationally invariant VBS. We would 
like to stress the fact that the validity of this picture is 
rooted in the existence of the Ni impurity. In fact, the 
physics of a homogeneous system of Cr atoms is closer 
to a gapless critical phase rather than to such a gapped 
state. 

General (to, n)-VBS states belong to a class of spin 
liquids which are known to possess an special hidden or- 
der that can be identified by a particular non-local or- 
der parameter called the string order parameter (SOP) 
[H, il, [13, m, m, [Hi . This order parameter has proved 
itself extremely successful in the task of characterizing 
diverse kinds of such states, both in the pure one dimen- 
sional cases and also in less trivial systems such as ladders 
[1^ , [131 . Wc shall see that this parameter also allows us 
to characterize the CrSNi ground state. The definition of 
the generalized string order parameter is as follows [2^ : 
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O,tr(0)= .lini {S^exp{i9Y,Sl)S;). (2) 



Ni Cri Cr2 Cr3 Cr4 Cr5 Cr6 Cr7 Cr8 

FIG. 3: Representation of two ways to measure the SOP in the 
CrSNi ring choosing adjacent starting sites. In the diagram 
at the top we break two bonds at the left end and one at the 
right. This corresponds to a (2, 1)-VBS. In the diagram at 
the bottom we have the same physical configuration but now 
we are breaking one bond at the left and two at the right end, 
which characterizes a (1,2)-VBS. Notice that both measures 
have been taken with an odd number of sites, otherwise the 
SOP is zero. 

The limit \ j — i\ going to infinite must be understood as 
comprising a region large enough so as to neglect finite- 
size efects. In our case the system itself is finite, however 
we will see that the results are still conclusive despite 
some corrections that have to be considered due to this 
fact. 

We hereby summarize the most relevant properties of 
the SOP for our purposes: 
i/ The SOP is symmetric with respect 9 = tt. 
ii/ The imaginary part of the SOP vanishes as we 
consider larger systems. 

iii/ Due to the fact that the definition of the SOP 
makes use of antisymmetric operators under spin flip, 
the number of sites considered in the measure must be 
odd in order to have an even number of these operators. 
Otherwise the SOP is zero. 

iv/ Given a generic (to,7i)-VBS state, the number of 
zeros of this operator in the interval 9 e [0, 27r) coincides 
with the number m. 

v/ Two measures of the SOP begining in adjacent sites 
will differ in the order of the numbers m and n. That is, 
if one measure gives a (to, n)-VBS state, the other will 
be a (n, m)-VBS (see fig. [3|. 

In order to show the existence of a VBS-picture of the 
CrSNi molecule as described in figO we have computed 
the generalized string order parameter in the ground 
state of the molecule for different ratios of the coupling 
constants J and J'. In virtue of property iv/ above we 
can characterize the VBS state counting the number of 
zeros of this operator. Moreover, we have some freedom 
at choosing the starting and end sites to measure the 
SOP (i.e, the numbers i and j in definition ([2])) and we 
will see that the results are consistent with property v/. 

In fig. m we have plotted the SOP computed in two dif- 
ferent blocks of the ring and adjacent starting sites (these 
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FIG. 4: a) SOP computed between the Cr2 and the Cre 
atoms, b) SOP computed between the Cri and the Cr^ 
atoms. We have fixed J — 16. Considering the interval 
£ [0, 2tt) and looking at the curves with high values of J' 
we see that the number of zeros in a) is two, corresponding 
to a (2, 1)-VBS. In b) the situation is more confusing due to 
stronger finite size effects that shifts the curves towards neg- 
ative values at — and 9 — tt. With this correction in mind 
the graph strongly suggests the existence of a (1, 2)- VBS. 



blocks are represented in fig. [3]). Remarkably, even in this 
small molecule, the imaginary part of the SOP vanishes 
in all our computations. 

In fig.S^ the curves with high values of J' have two local 
minima in the interval 9 € [0, 27r) whose value is compat- 
ible with zero considering that the system is finite and 
the block used to measure the SOP comprises only some 
few spins. This result is consistent with a (2, 1)-VBS. 
The shape of the SOP is also typical of these states, with 
two noticeable maxima placed approximately at 9 = n/A 
and 9 ~ 37r/4. Accordingly, the computations shown in 
fig.lHa shall be consistent with a (1, 2)-VBS as explained 
above. We see that again the shape of the SOP for the 
strong coupling curves shares the main features of these 
states, that is, one substantial maximum placed precisely 
at 6* = tt/2. There is however a major difference at = 



5 



and 9 = IT, where the SOP does not vanish but has a sig- 
nificant negative value. At these points the exponentials 
in the definition of the SOP do not play any role and 
hence it is a usual spin-spin correlator between the two 
Cr atoms at the end of the block. Since our ring is finite 
and closed it is quite natural to think that the Cri and 
Crr spins are highly correlated via the spin frustration in 
the Ni impurity. Therefore, such a finite value at = 
only shows the compromise between the bulk (1, 2)-VBS 
picture of the whole system and the local physics of those 
spins. 

In summary, the string order parameter reveals that 
the ground state of the CrSNi in the strong impurity cou- 
pling limit is consistent with a impurity-mediated mech- 
anism of spin frustration, the result of which is the VBS- 
pattern shown in fig.[H Since we are dealing with a finite 
system the possibilities to drastically change the nature 
of the ground state are limited, that is, the possibility of a 
quantum phase transition is excluded and only a crossing 
with another energy level can produce this effect. There- 
fore, the intermediate J'/ J regime can be considered as 
some deformation of the strong coupling limit. As we 
decrease J' and the Stot = 1 state crosses the Stot = 0, 
the VBS picture breaks down and the measures of the 
SOP are not meaningful in the sense that the properties 
of this operator in such a state are not well defined. 



havior of the CrSNi molecule by means of computing the 
time evolution of some important spin correlators: the 
spin autocorrelation of the impurity Ni atom with itself 
and the spin correlation between the Ni atoms and each 
Cr along the ring. These correlators correspond to the 
vacuum expectation value of the time-evolved spin oper- 
ators SNi(i) and Scii(i) projected onto the spin operator 
of the Ni impurity at t = 0, SNi(0). This is a way to 
dynamically probe [H, [2^, [30, [Ml, [111 the spin structure 
in the ground state |^o) of the ring molecule. In fact, 
we shall consider the square modulus of those correlators 
and interpret them as time-evolution probabilities. That 
is, we shall consider the following correlators in order to 
construct a figure of merit: 

CNi(i) := (V'olSmW • SNi(0)|?/'o), (3) 

and 

Cc,,W :=(Vo|Scr.W-SNi(0)|Vo), (4) 

Notice that in a rotationally invariant singlet ground 
state \'ipo) the correlations in the x,y and z axis have 
the same value and thus, the expresions above can be 
written as: 

CN^{t) - HMSk{t)sum^), (5) 



IV. FRUSTRATION AND DYNAMICS IN 
CRSNI 

As we have already mentioned, the CrSNi ring molecule 
is frustrated in the sense that the minimum energy of the 
system can not be obtained minimizing separatedly each 
of the two body terms of the Hamiltonian ([1]). Another 
way to see it is by resorting to the classical limit where 
each spin is pictured as a classical vector. Once we set 
the value of one of those classical spins, then we can fix 
one by one the rest of the spins in order to minimize the 
local interactions, but in the end there will be one spin 
for which the local interaction with both of its neighbors 
can not be minimized at the same time. 

Typically, systems where frustration exists come along 
with a rich and very often not so well-known physics. 
Roughly speaking, we can say that frustration in general 
increases the complexity of those systems, both in the 
physics they exhibit as well as in the way to approach 
them. In particular, there is not a well defined way to 
measure the amount and localization of frustration. An 
attempt to quantify these effects in CrSNi can be done 
attending to the structural changes of the gound state 
as we vary the couplings. That is, by inspection of the 
way in which spins in the ring couple to form the final 
state. This procedure has a connection with experimen- 
tal techniques where the Lande factors of the ring can 
be measured. However this procedure is not suitable to 
study a rotationally invariant singlet ground state where 
the spin is zero. In this section we will study the be- 



CcnW = 3(^o|5a,,(i)5Ni(0)|^o), (6) 

Since the proportionality factor does not provide any ad- 
ditional information we will discard it from now on and 
will consider the bare z-axis projection correlators. The 
time dependency of the operators is given by the usual 
Heisenberg picture: 

0{t) = e'"'O{0)e-'"\ (7) 

The idea behind this figure of merit to measure the 
dynamical correlations between spins is similar to the 
static correlator used to measure spin correlations in 
space separated sites i and j of the ring (^/lolScri (*) ■ 
SNiO)|V'o)- This static correlator measure spatial cor- 
relations, while our purpose is to measure time-evolved 
correlations which will probe not only the ground state 
physics but also the excited states physics. 

We next explain briefiy the numerical method used to 
evaluate these correlators. After that we shall show and 
discuss the results. 



A. Numerical Method 

The Hamiltonian ([T]) is SU(2) rotational invariant. 
That is, it commutes both with the total spin and the 
z-axis projection of the total spin. In table |T| we show 
the dimensions of the subspaces corresponding to the 
conserved quantum numbers of these operators. These 
sizes are in the limit to perform exact diagonalization 
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Dimension 


Dimension 





1000 


23548 


1 


2764 


22548 


2 


3905 


19784 


3 


4256 


15879 


4 


3900 


11623 


5 


3095 


7723 


6 


2150 


4628 


7 


1308 


2478 


8 


692 


1170 


9 


314 


478 


10 


119 


164 


11 


36 


45 


12 


8 


9 


13 


1 


1 


Total: 


196608 


196608 



TABLE I: Dimension of each subspace with well defined 
quantum numbers out of the total Hilbert space of a CrSNi 
ring. The second column corresponds to the subspaces with 
well defined total spin. The third one are the sectors with 
well defined value of the z-axis projection of the total spin. 
In this last case for each value in the first column we must 
consider the positive and negative cases. 

but lie however within the domain reachable for a Lanc- 
zos method. 

To do our computations we will make use of an adapted 
version of the Lanczos al gori thm specific to compute real- 
time dynamics [33l [3^ . l35j | . In this framework it can 
be shown that the correlators written above can be ex- 
pressed as: 

M 

3=0 («) 

where M stands for the dimension of the Krylov space 
]C{H,qo,M) such that K{H,qo,M) = /C(iJ,(?o,M + 1), 
with go := Sf^iH'o)- That is, M is the dimension of 
the largest invariant subspace generated by succesive 
aplications of the Hamiltonian H upon the seed vector 
go = SNi\'^o)- The vectors j^/'j) are the aproximated 
eigenvectors computed in this Krylov subspace, are 
the energies of these eigenvectors and Eq stands for the 
energy of the ground state. 

The number M is typically much lower than the total 
dimension of the Hilbert space, but still high to numer- 
ically compute a complete basis of the Krylov subspace. 
Therefore, the approximation in this method resides in 
the fact that we will substitute the dimension M with 
a lower number of vectors that still serve as a complete 
basis for these correlators. 

In order to obtain the most accurate results and repre- 
sentations of the eigenvectors of the Hamiltonian we 
have not used the same Krylov space to compute them 
all. Instead we have performed a Lanczos iteration to 
find the ground state. After that, the Lanczos iteration 



is restarted with the previously found eigenvectors pro- 
jected out of the subspace to find the next excited state. 
And so on and so forth, until we have computed the de- 
sired number of eigenvectors. In particular, to compute 
the correlators described before we have used 400 eigen- 
vectors of the Hamiltonian. As regards of the tolerance 
in the eigenvalues, we have set it to 10""'^"' allowing a 
maximum dimension of each Krylov space of 350 vectors. 
Were we in an exact situation, these vectors should be 
normalized to one and be orthogonal among themselves. 
Let us call V the matrix whose columns are these eigen- 
vectors. We have checked that we obtain the following 
accuracy, 

\\Wv-i\\^iir\ (9) 

which can be considered a low value for such a large 
number of eigenvectors. Moreover, as another check of 
the accuracy of the eigenvectors we have computed the 
total spin of each one of them and we have obtained inte- 
gral values up to precisions of 10^^ in the vast majority 
of them. As for the invariance of the Hamiltonian un- 
der reflection with respect the Ni atom, we have checked 
that symmetric one and two body correlators evaluated 
on every eigenvector are the same up to the fourth or 
fifth decimal digit. 



J' 


Cm[Q) 


Error (%) 


2 


0.66666 


0.00002 


10 


0.66666 


0.0009 


20 


0.66653 


0.02 


30 


0.66512 


0.2 


40 


0.65660 


2 


50 


0.57469 


14 



TABLE II: Relative error between the Ni self-correlation 
at t=0 using Q with M=400 eigenstates and the value 
{il)o\S^iS^i\ipo) in the ground state, which is equal to 2/3 . 

In order to check how complete is our set of eigenvec- 
tors, we have compared the value at < = computed 
using ^ and (g]) with the values obtained measuring the 
correlators in the ground state without the projectors in 
between. We have observed that the agreement is exce- 
lent for low values of J' while it goes worse for higher 
values of the coupling constant. Table [IT] shows these 
values and the relative error. 

These accuracy checks confirm that our numerical re- 
sults are good enough so as to interpret them on physical 
grounds with respect to the spin frustracion effects de- 
scribed in previous sections. 



B. Results 

In the previous sections we have proposed and checked 
with the proper order parameters a static picture of the 
ground state of the CrSNi ring. In the limit when the Ni 
is weakly coupled to the Gr bulk {J' / J <Si 1) the system 
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FIG. 5; Modulus of the spin correlation of the impurity Ni 
atom with itself and with the rest of the Cr atoms. The 
graphs have been done with a value J = 16 K and: a) J' — 2 
K, b) J' = 10 K, c) J' = 20 K. For convention the sign of the 
correlators has been chosen to coincide with the sign of the 
correlation at t=0, which is a real number. 



FIG. 6; Modulus of the spin correlation of the impurity Ni 
atom with itself and with the rest of the Cr atoms. The 
graphs have been done with a value J = 16 K and: a) J' = 30 
K, b) J' = 40 K, c) J' = 50 K. For convention the sign of the 
correlators has been chosen to coincide with the sign of the 
correlation a,t t — 0, which is a real number. 



is accurately described by an isolated Ni atom and an 
open Cr8 chain. On the other hand, when the impurity 
is strongly coupled {J' /J 3> 1) the ring posess a dom- 
inant contribution in the form of a (2, 1)-VBS ground 
state with some local correlations around the Ni atom 
due to the finite size of the sample. From the point of 
view of frustration, in this limit the impurity acquires 
strong antiferromagnetic compromises with both of its 
neighbors that can not satisfy simultaneously. Frustra- 
tion is known to impose complex constraints that can 
destabilize, deform and even produce new states of mat- 
ter. On the other hand, the dynamics of each spin of the 
system is highly influenced by these constraints. In the 
following paragraphs we will see that the self correlation 
of the Ni impurity and the rest of correlators with the 
Cr atoms allows us to naturally establish a relation with 
the amount of frustration. 

In figs. [5] and [S] we have plotted the correlators (O 
and (2]) for a fixed value of the constant J = 16 K and 



different values of J'. The CrSNi ring is invariant un- 
der reflection respect the impurity site and therefore we 
will only provide the correlators with the Ni itself and 
the Cr atoms numbered from 1 to 4 (with the notation 
of fig. [2|) , the correlations with the Cr atoms numbered 
from 5 to 8 are the same as their symmetric counterparts. 
We want to provide these magnitudes with the meaning 
of a time-evolved probability and hence we will consider 
only their modulus. It is worth noticing that the time 
correlators mentioned before, at i = are real numbers 
whereas for arbitrary values of t they are complex num- 
bers. For convention, in the graphs where we plot the 
modulus of these correlators we will provide them with 
the same sign of their real value at < = to make explicit 
the ferromagnetic or antiferromagnetic nature that they 
possess in the static t = ground state. 

In fact, we can observe in these graphs that the fin- 
gerprint of an antiferromagnetic order is present in the 
initial t = ground state and shows up in the alternation 
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FIG. 7: The number of oscillations has been obtained count- 
ing the number of minima in a wide time interval and divid- 
ing by the total time. The Ni, Cr\ and Cr2 have increas- 
ing frecuencies with J' which corresponds to more frustrated 
regimes. In particular for the Ni spin the relation of these 
variables is linear up to a high precision. In the case of the 
Cra and Cr^ spins the frecuency is hardly affected by J'. 



of the signs of the correlators. Notice also a signal of frus- 
tration in the fact that this alternation fails in the C7'4 
and its symmetric counterpart Crs (due to the reflection 
invariance of the ring their value is equal with the same 
sign), where the correlations reveal that both spins are 
oriented in the same direction respect the spin of the Ni. 
Remarkably this ferromagnetic defect is a consecuence 
only of the reflection invariance of the Hamiltonian. 

From these graphs we can also infer that the average 
correlation of each spin with the impurity is little sensi- 
tive to the strength of the coupling constant J', although 
the amplitude of the deviations with respect to this av- 
erage value increases with it. 

The most important observation is that the dynamics 
of the correlators exhibit a non trivial sort of periodicity. 
That is, from the shape of the curves it seems that there 
exists many modulating components but a dominant pat- 
tern of oscillations is apparent. Moreover, the frecuency 
of this pattern clearly increases with increasing values of 
the impurity coupling J', i.e. as we move towards more 
frustrated regimes, but not in the same way for all the 
spins. We have captured in figl7| the frecuency of the 
dominant oscillatory pattern of each correlation. The 
graph highlights two different tendencies of the correla- 
tions depending on the considered spins: the frecuency 
of the Ni self correlation as well as the correlations of 
the Ni spin with the Cri and Cr2 spins increase with 
J'. Moreover in the case of the Ni self correlation the 
relation of these two variables is linear with a surprising 
accuracy. On the other hand the Cr-^ and Cr^ are less 
affected by the impurity spin and the frecuency remains 
almost constant in the wide range considered. 

In fig. |5] we have plotted the same real time correla- 
tions with an election of the coupling constants J = 16 
K and J' = 1 K such that the ground state is a triplet. 
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FIG. 8: Modulus of the spin correlation of the impurity Ni 
atom with itself and with the rest of the Cr atoms. The values 
of the coupling constants are J = 16 K and J' = 1 K which 
corresponds to the region where the ground state is a triplet. 
For convention the sign of the correlators has been chosen to 
coincide with the sign of the correlation at t = 0, which is a 
real number. 



In this case the ground state is not rotationally invariant 
and the z-axis projection of the correlators is not pro- 
portional to the scalar correlators ([3]) and ^ . For our 
purposes however, these magnitudes suffice to realize the 
different nature of both the singlet and triplet ground 
states: first of all is that not one but two dominant pat- 
terns of oscillation are well distinguishable in the triplet 
regime. Secondly and also a major difference is that the 
scale in this regime is some orders of magnitude smaller 
than the singlet case. 

The existence of this oscillatory behavior seems natu- 
ral in a frustrated system where there does not exist a 
natural equilibrium position for each spin or where the 
resulting equilibrium configuration may result unstable. 
In a system composed by classical spins these oscilla- 
tions can be interpreted as the necessary movements of 
each spin to satisfy the frustrated interactions, becoming 
faster as we blur the concrete equilibrium positions with 
the frustrating interactions. 

The results in fig. [7] point towards a regime where 
the frustration introduced by the Ni impurity has strong 
local dynamical effects in the nearest and next nearest Cr 
neighbors while the rest of the spins perceive the impurity 
screened by this closer shell of atoms and therefore their 
dynamics is little affected by it. These results also show 
that the correlators proposed to study the frustration 
of the system indeed have the behavior expected for a 
suitable estimator in order to measure the intuitive idea 
we have about the amount of frustration in a certain 
system. 



V. CONCLUSIONS 

In recent years, considerable efforts have been devoted 
to synthesizing and investigating magnetic systems of 
nano scale dimension that comprise a controllable num- 
ber of transition metal ions. Highly symmetrical clus- 
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ters of almost planar ring shape are among such topi- 
cal molecular nanomagnets. In particular, the bimetallic 
ring molecule CrSNi is the first antiferromagnetic ring 
with an odd number of spins. Thus, it is a remarkable 
quantum system to test fundamental magnetic proper- 
ties, and in particular the spin frustration effects. 

In this work, we have studied the CrSNi frustrated ring 
in the regime where the ground state is a singlet. That 
is, with a fixed value J = 16 K this region corresponds 
to J' > 1.5 K. In this regard, the experimental charac- 
terization of a CrSNi molecule places the real strengths 
present in the real system close to J = 16 K and J' = 70 
K, well within the singlet region. 

As we let the interaction strength of the Ni impurity 
to be stronger than that between the Cr atoms, the ring 
stabilizes in a ground state with the quantum proper- 
tics of a dimerized VBS. The picture that explains this 
behavior in terms of the possible bonds between neigh- 
bor particles comes clear from fig. [21 Such a VBS state 
constitutes an example of a spin liquid with an intrinsic 
order that can be measured by means of some particular 
non-local order parameters. In fig|3]wc show the compu- 
tations of this order parameter on the ground state and 
its behavior supports neatly the VBS picture. In this re- 
gard, some finite size effects can be observed in the order 
parameter that reveal a competition between the physics 
in the bulk of the ring and the strong effects, possibly 
mediated by the system frustration, that ocurrs in the 
vicinity of the Ni atom. 



In the second section of this paper we have studied the 
role of the frustration in such a VBS state by means of 
computing the real time evolution of the spin correlators 
between the atoms in the ring. In particular, we have 
found that the amount of frustration can be related to 
the frecuency in the oscillatory behavior of this correla- 
tors. This relation can be naturally established from the 
observation that the oscillations in the system become 
faster as we move to the more frustrated regime J' ^ J. 
Such an oscillatory behavior is natural in a system where 
no natural equilibrium is allowed due to the frustration. 
However, the spin correlators reveal that the atoms that 
are most affected by this frustration are the Ni impurity 
itself and those Cr atoms that are closer to it, that is Cri 
and Cr2, while the effect of the impurity strenght seems 
to be less influent in the Cra and Cr4 atoms. 

We believe that the methods and numerical techniques 
used in this work are versatile enough and can be ex- 
tended to a variety of other nanomolecular magnetic com- 
pounds. 
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